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ABSTRACT 

The size distribution of the stability region around the Lagrangian point Li is in- 
vestigated in the elliptic restricted three-body problem as the function of the mass 
parameter and the orbital eccentricity of the primaries. It is shown that there are 
minimum zones in the size distribution of the stability regions, and these zones are 
connected with secondary resonances between the frequencies of librational motions 
around _L 4 . The results can be applied to hypothetical Trojan planets for predict- 
ing values of the mass parameter and the eccentricity for which such objects can be 
expected or their existence is less probable. 
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1 INTRODUCTION 

The existence of Trojan planets in extrasolar planetary sys- 
tems (EPS) is an interesting possibility. Some related prob- 
lems have been recently studied in several papers. In the 
Solar system, more than 2000 Trojan asteroids have already 
been discovered. They move close to the orbit of Jupiter, 
• near the stable Lagrangian points L4 and L5, and are in 
1:1 mean motion resonance with the giant planet. Mars 
and Neptun also have co-orbital asteroidal companions, and 
there are co-orbital moons around Saturn as well. It can be 
expected that Trojan-type objects also exist in EPSs. 

So far there have been discussions on Tro jan exoplan- 
ets. According to lLaughlin fc Chambers! (|2002l ). two planets 
around a solar-mass star, with masses comparable to the 
mass of Jupiter or Saturn, can perform stable tadpole-type 
librations around the Lagrangian points. Pairs of Saturn- 
mass planets can also execute horseshoe orbits around a 
solar- mass star, but th is is not possible for Jupiter-mass 
pairs. iNauenberd |2002l ) determined numerically the non- 
linear stability domain of the triangular Lagrangian solu- 
tions in the general three-body problem as a function of the 
eccentricity of the orbits and the Routh's mass parameter. 
This study indicates that there is a wide range of Jupiter-size 
planetary masses and eccentricities for which such solutions 
could exist in EPSs. 

As to the formation of Trojan planets, several 
mechanisms have been pr oposed, including ac c retion 
from protoplanetary di sk l|Laughlin fc Chambers! 120021: 

IChiang fc Lithw ick 2005]) and capture into 1:1 resonance 
due t o rapid mass growth of the planet l|Chiang fc Lithwick! 

120051 ). or convergent migration of multiple protoplanets 
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Thommes! I2005I : | Cresswell fc Nelson! 120061 ). iBeauge et al.l 
20071 ) studied the formation of terrestrial-like Trojan plan- 



ets from a swarm of planetesimals around the Lagrangian 
point L4 of a giant planet. For the form ation of the Tro- 
j an as teroids of Jupiter and Neptune, see iMorbidelli et al.l 
l|2005l ). 

Trojan exoplanets could be detected in several ways, 
like radial velocity (RV), astrometric, or transit photom- 
etry measurements. lLaughlin fc Chambers! |2002i ) pointed 
out that a pair of comparable mass planets, both in 
tadpole and horseshoe orbits, would induce a character- 
istic pattern in the RV component of the central star 
that could be detected. Reanalysing the RV measurements 
of the systems HD 128311 and HD 82943, where two 
Jupiter-like planets are beleiv ed to be in a 2:1 resonance, 
IGozdziewski fc Konackil l|2006h demonstrated that it is pos- 
sible to explain the variations of the RV of the hosting stars 
by a pair of Trojan planets in 1:1 resonance. According to 
IGozdziewski fc Kon acki (2006), this is also possible f or the 
syst em HD 73526. However, in t his ca se iTinnev et ail (|2006l ) 
and ISandor. Kiev fc Klagyivikl (|2007l ) found only a 2:1 reso- 
nance. |Foid_&_CJ^udj (|2006l ) developed a method for detect- 
ing Trojan companions to extrasolar planets by combining 
RV and photometric observations of transiting exoplanets. 
This method offers the potential to detect even terrestrial- 
mass Trojan planets using ground-based observations. 

Most of the known exoplanets are gaseous giants or- 
biting quite close to their hosting stars. The search for 
terrestrial-like planets is the ambitious aim of several ongo- 
ing and future ground-based and space projects (COROT 
has been in orbit since 2006 December, Kepler will be 
launched in 2008, others are planned for the next decade). 
It is a question of considerable interest, whether Earth-like 
planets exist in the habitable zone (HZ) of some exoplane- 
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tary systems. If there is a giant planet in the HZ of a sys- 
tem, the exi stence of another planet th ere is unprobable. 
However, as iMenou fc Tabachnikl (|2003l ) noted, terrestrial 
planets could exist near the stable Lagrangian points L4 or 
L5 of a giant planet movin g in the HZ. 

lErdi fc Sand"oil (|2005h studied this possibility in de- 
tail, investigating five exoplanetary systems (HD 17051, HD 
28185, HD 108874, HD 177830, and HD 27442) in which 
the only known giant planet moves in the HZ. By us- 
ing the model of the elliptic restricted three-body problem 
(ERTBP), they numerically determined the region around 
Z/4 of each system where stable tadpole-type librations are 
possible. Four other systems (HD 150706, HD 114783, HD 
20367, and HD 23079) were also studied in which the orbit 
of the giant planet is partly outside the HZ due to its large 
eccentricity. It has been shown that in all studied systems 
there is an extended stability region around L4, whose size 
depends on the mass ratio of the planet and the star and 
the orbital eccentricity of the giant planet. It is possible that 
Trojan exoplanets of negl igible mass exist in these systems. 
Othe r dynamical studies l)Dvorak et al.ll2004l ; ISchwarz et al.l 
120051 ) also confirmed the possibility that small Trojan plan- 
ets can exist in stable orbits with moderate eccentricities in 
thes e systems. 

ISchwarz et all l|2007f ) reconsidered the system HD 
108874, where a second giant planet was discovered 
l|Vogt et al.ll2005l ) outside the HZ, perturbing strongly the 
inner planet. They showed by numerical integrations for 10 7 
yr that for most combinations of the orbital parameters, 
motions around L4 are unstable. However, in some cases, 
terrestrial-like Trojan planets can exist in a small domain of 
stability. 

Hypothetical Tro jan plan e ts in other multiple EPSs 
were also considered. I Ji et alj l|2005h showed that Trojan 
planets with 0.1 - 10 Earth-mass and with low eccentrici- 
ties and inclinations could stay for 10 7 yr at the triangular 
equilibrium points of the tw o massive planets of 47 UMa. 
According to I Ji et al.l (|2007i ). Trojan planets with 0.01 - 1 
Earth-mass can also exist at least for 10 6 yr in the triple 
planetary system HD 69830 in 1:1 mean motion resonance 
with the outer massive planet d. 

In this paper, we investigate how the size of the sta- 
bility region around L4 depends on the parameters of the 
problem. Since in more than 90% of the EPSs discovered so 
far only one (giant) planet is known, we use a model of the 
three-body problem corresponding to a gravitational system 
of a star, a giant planet and a Trojan planet. We study cases 
where the Trojan planet has negligible mass and the orbit 
of the giant planet is elliptic (ERTBP), and also the Trojan 
planet can have finite mass and then we assume that the 
orbit of the giant planet is a perturbed ellipse. We deter- 
mine the size of the stability region around La depending 
on the mass of the giant planet (with respect to the mass 
of the star) or the masses of both planets and on the (ini- 
tial) orbital eccentricity of the giant planet. For the sake of 
simplicity, here we only study coplanar orbits of the planets. 



2 MAPS OF DYNAMICAL STABILITY 

To determine the size of the stability region around L4, we 
have to establish the dynamical character of orbits in the 
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Figure 1. The stability regions around L4 and L5 for e' = 
and fi = 0.001. There can be temporary stable horseshoe orbits 
outside the main stability regions where orbits have librational 
character. 



neighbourhood of this point. In the model of the ERTBP, for 
this we used the Lyapunov characte ristic indicators (LCI) 
l|Benettin et al.lll980l ; lFroeschlelll984h . computed for 10 4 pe- 
riods of the giant planet. During the numerical integration, 
we also checked the maximum eccentricity of the Trojan 
plane t as another tool o f dynamical characterization of or- 
bits l|Dvorak et all koOS). We used dimensionless quantities 
with properly selected units. The sum of the masses of the 
star (mo) and the giant planet (m') was taken as unit mass. 
Then the mass of the star and the giant planet is 1 — p and 
fx, respectively, where fi = m' /(mo + m') is the mass param- 
eter. The semi-major axis a of the orbit of the giant planet 
was taken as unit distance. 

We computed maps of dynamical stability depending 
on two parameters, fi and the eccentricity e' of the orbit 
of the giant planet, that were changed in specific intervals. 
We changed /1 between 5 x 10~ 4 and 1.5 x 10~ 2 in steps 
of A/i = 5 x 10 -4 . (With a one solar-mass star this corre- 
sponds to a giant planet with about 0.5 - 15 Jupiter-mass.) 
The eccentricity e' was changed between and 0.6 in steps 
of Ae' = 0.025. For each possible pair of /i and e', we com- 
puted a stability map showing the stable and unstable re- 
gions aroud L4 in the parameter plane (a, A — A'), where a 
is the semi-major axis of the orbit of the Trojan planet and 
A — A' is the synodic longitude, A and A' being the mean 
orbital longitude of the Trojan and the giant planet, respec- 
tively. 

To obtain a stability map, we changed a between 0.9 
and 1.1 in steps of Aa — 0.0025, and A between 0° and 180° 
in steps of 1°, taking the direction of the pericenter of the 
giant planet as starting point. The other initial conditions 
were: A' = 0°, e = 0, where e is the eccentricity of the orbit 
of the Trojan planet. At each point of the (a, A — A') plane, 
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with the given initial data we computed the value of the 
LCI for the resulting orbit and represented its logarithm on 
a colour-coded scale. In what follows we discuss the main 
characteristics of the computed maps. Some representatives 
of them are shown in Figs 1-3. Light regions correspond 
to stable orbits (at least for the investigated time interval). 
The dark background corresponds to escape or collision or- 
bits with the giant planet. Computations of the maximum 
eccentricity show that in such cases the eccentricity of the 
Trojan planet grows above 0.3. 

Fig. 1 shows the stability regions around L4 and L5 
for e! — (circular restricted three-body problem) and n = 
0.001. To compute this figure, we additionally changed A 
between 180° and 360°. In the two main stability regions 
orbits have librational character. Outside the main stability 
regions there can be temporary stable horseshoe orbits. 

Increasing the mass of the giant planet, the stability 
region becomes shorter in the synodic longitude and wider in 
the semi-major axis. In the bottom panel of Fig. 2, obtained 
for e' = 0, fi = 0.002, a ring structure appears near the 
edge of the stability region, corresponding to a high order 
resonance between the short and long period components of 
librational motions around L4. Outside the main stability 
region, a chain of small islands can be seen, originating from 
the disruption of a former ring, due to the increase of the 
mass and increasing perturbations of the giant planet. 

One expects that with the increase of the mass of the gi- 
ant planet, the extension of the stability region shrink due to 
larger perturbations. However, the comparison of the middle 
and top panels of Fig. 2 shows the opposite. In the former 
case (e' = 0, fj. = 0.006), the stability region is smaller than 
in the latter for larger fj, (e! = 0, n = 0.007). This is one 
result of our computations: with the increase of the mass 
of the giant planet, the size of the stability region around 
L4 does not decrease uniformly, smaller and larger regions 
alternate along the overall decreasing tendency. 

This phenomenon can be also seen in the ERTBP. Fig. 3 
shows for e' = 0.1 the stability region for increasing values of 
[l. Its size is decreasing as \x takes increasing values (bottom 
and middle panels) , and it is larger again for a higher value 
of fl . 

Studying the computed stability maps (not shown 
here), we can arrive at a similar conclusion: when we fix 
[i and change e', with the increase of the eccentricity the 
decrease of the size of the stability region is not uniform. 

We also studied cases when the Trojan planet 
had non-negligible mass. Because of the larger complex- 
ity of the problem, in such cases we used the effi- 
cient method of the relative Lyapunov indicators (RLI ) 
|Sandor. Erdi fc Efthvmiopoulosl l2000l : ISandor et all I2004T I 
to determine the dynamical character of orbits around L4. 
We computed stability maps on the (a, A — A') plane by 
considering several values for the mass of the Trojan planet, 
corresponding to m = 1, 10, and 100 Earth-mass (the) and 1 
Jupiter-mass (m,j), and for the mass of the giant planet, cor- 
responding to 1, 2, ... , 7mj, assuming that the planets move 
around a solar-mass star. For the initial eccentricity of the 
giant planet, values were taken from the interval e' = 0—0.30 
in steps of 0.05. The other initial conditions were as before. 

The results of our computations show that the extension 
of the stability region do not change much under the effect 
of the mass of the Trojan planet. Even for m — 100 ms, the 
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Figure 2. The stability region around L4 for e' = and fi = 
0.002 (bottom panel), = 0.006 (middle panel ), and fi = 0.007 
(top panel). Note the ring structure near the edge of the stability 
region in the bottom panel and the chain of small islands in the 
bottom and middle panels. 



stability regions look like in Fig. 3 (for the corresponding 
values of m' and initial e'). 

We determined the stability regions around L4 for the 
posssible combinations of the values of m, m' and e' taken 
from the previously mentioned sets. The results show that 
for a given pair of m' and e', the size of the stability region 
does not change much with the increase of m. The changes 
are larger when m is fixed, and either m' or e' is changed 
while the other is kept constant. The non-uniform decrease 
of the size of the stability region is observed in such cases as 
well. The computations confirm the existence of a stability 
region around L4 even for m = mj , when the giant planet is 
of several Jupiter-masses and its orbit is very eccentric. As 
an example, Fig. 4 shows the stability region for m = lmj, 
e' = 0.3, and m! — lmj. Increasing m' at this value of m 
and e', the size of the stability region decreases reaching a 
minimum at m' = 5mj, after which it is larger again for 
m' — 6mj, as the computations show. (Note, however, that 
in this case the stepsize is only Am' = lmj.) 

From the numerical investigations of this section, two 
main conclusions can be drawn: (i) the size of the stability 
region around L4 depends essentially only on the mass pa- 
rameter n and the eccentricity e'; (ii) the variation of the 
size of the stability region is not uniform either with chang- 
ing fi or e'. To find an explanation for the second property 
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Figure 3. The stability region around L4 for e' = 0.1 and fi = 
0.002 (bottom panel), fi = 0.007 (middle panel), and fj, = 0.009 
(top panel). 




Figure 4. The stability region around L4 in the planar three- 
body problem for m = lmj, e' = 0.3, ml = lmj. 



needs a more thorough study, since though we considered 
many values of \x and e , however, the applied resolution is 
not enough to get a satisfactory description of how the size 
of the stability region depends on these parameters. In the 
next section, therefore, we determine the distribution of the 
size of the stability region in a simplified model with high 
resolution as the function of /1 and e'. This will enable us 
to study the fine structure of the size distribution of the 
stability region. 



3 SIZE DISTRIBUTION OF THE STABILITY 
REGION AROUND L 4 

The extension of the stability region around the Lagrangian 
points L4 and L5 is an old problem of celestial mechanics 
that was studied in different contexts by many authors. Most 
of the papers refer to the problem of the stability of Jupiter'; 
Troja n asteroids (see e.g. Levison. S hoemaker & Shoemaker 



19971; iGiorgilli fc Skokosl [l997l; lEfthvmiopoulos fc Sandoi 



20051 ; fRobutcl & Gabcrn 20061 ). 



In the model of the ERTBP, the stability of the point 
L4 depends on the mass parameter [i and the orbital ec- 
centricity of the primaries which from now on we denote 
by e. The stability of orbits around L4 also depends on 
these paramete rs. Non-linear st ability of orbits around L4 
was studied by iGvorgvevI (|l985l ) for the mass parameter of 
the Ea rth-Moon system and for a few values of the eccen- 
tricity. lLohinger fc Dvorak! <| 19931 ) investigated the size of 
the stability region around L4 by computing the number of 
stable orbits for a wide range of the values of /i and e. The 
most detailed survey on the size distribution of the stabil- 
ity region around L4 depending on /x and e was given by 
lErdi fc Sandoi] (|2005l ). 

Here we reconsider the problem by using a much higher 
resolution in fi and e than in previous works. Fig. 5(a) shows 
the size distribution of the stability region around L4 in the 
ERTBP depending on the mass parameter /1 and eccentricity 
e. Note, that here and from now on /i has slightly different 
meaning than in Section 2. In order that we could compare 
our results with earlier works, we changed the system so 
that now the larger primary (the star) has a mass 1 and 
the smaller primary (the giant planet) has mass /1 (thus 
their total mass is 1 + fi, instead of 1). The test planet (the 
Trojan planet) has negligible mass. 

Fig. 5(a) was obtained as follows. For a given pair of 
the values of and e, we placed the test planet in the point 
L4, and moving it outwards along the line connecting L4 
and the larger primary, we checked by numerical integra- 
tion for 5 x 10 4 periods of the primaries the stability of the 
orbits, originating from consequtive points of this line with 
zero initial velocity with respect to the primaries. We de- 
fined stability in the simplest way: an orbit around L4 was 
considered stable if it did not cross the line of the primaries. 
The reason for this definition was that when a crossing hap- 
pens the orbit loses its librational character and we were 
interested in tadpole-type librations as possible orbits for 
Trojan planets. For any pair of fj, and e, we determined the 
largest distance from L4, denoted by e, along the line of L4 
and the larger primary, until the orbit of the test planet 
starting from points of this line remains stable in the above 
sense. Certainly, the time interval in stability investigations 
is crucial. We know that the exte nsion of the stabilit y re- 
gion around L4 shrink with time l|Levison et al.|[l997l ). We 
took 5 x 10 4 periods of the primaries for our stability inves- 
tigations as a compromise. The general features of the size 
distribution of the stability region around L4 as the function 
of fi and e appear reliably during this time. On the other 
hand, the high resolution computations could be performed 
in a reasonable time. 

Changing /1 and e on a fine grid (0.0001 ^ /1 ^ 0.05, 
Afi = 0.0001, < e < 0.6, Ae = 0.001), we determined for 
each pair of \i and e the limit value e (in the fraction of the 




initial distance of the primaries) corresponding to the largest 
stable librational orbit. The distance e can be thought as the 
approximate half thickness of the stability region described 
by the largest stable librational orbit. 

Fig. 5(a) shows the distribution of e as the function of 
H and e, that is the size distribution of the stability region 
around L4 as the function of the parameters. For the sake of 
better visualization, instead of e the values of l/log(e) are 
plotted in Fig. 5(a). The light region above the V-shaped 
boundary curve and on the right of the figure corresponds 
to instability, even the point L4 is unstable for such values 
of fi and e. This boundary curve is well known from stability 
investigatio ns of the point L4 in the ERTBP (|Danbvlll964l; 
lMeirelll98lD . The point L 4 is stable and libration around L4 
is possible only below this curve. The novelty of Fig. 5(a) is, 
that in the region where L4 can be stable, it indicates the size 
of the stable librational region around L4. According to the 
colour code, darker domains correspond to more extended 
stability regions. 

As can be seen from Fig. 5(a), the size distribution of 
the stability region has a very complex structure. The size 
of the stability region is the largest when both fi and e are 
small (/1 < 0.01, e < 0.1). The fine structure of the figure 
confirms our previous finding that the size of the stability 
region changes non-uniformly either fixing e and varying /i, 
or vice versa. It can be seen that there is also an extended 
stability region for small values of the eccentricity (e < 0.1) 
between ^ = 0.014 - 0.020. 

The dominant features of Fig. 5(a) are the dark and 
light stripes, forming up finger-like structures on the left of 
the figure. One might suspect that the light stripes, indicat- 
ing less extended stability regions, are related to resonances 
between the frequencies of librational motions. In the next 
section, we consider this possibility and see how resonances 
of librational motions around L4 in the ERTBP can be de- 
termined. 

We also computed the size distribution of the stability 
region around L4 by giving masses to the Trojan planet 
and we found that the structure shown in Fig. 5(a) remains 
even for a Trojan planet as large as 1 Jupiter-mass. The only 
difference is that with the increase of the mass of the Trojan 
planet, the stripes of the figure are slightly shifted toward 
the direction of lower mass parameters. Also, there is only 



a slight difference in the size distributions corresponding to 
pericentric or apocentric (as in Fig. 5a) initial positions of 
the primaries. Thus the structure of Fig. 5(a) seems to be 
quite general. 



4 RESONANCES OF LIBRATIONAL 
MOTIONS 

In the circular restricted three-body problem, motions 
around L4 has two frequencies, n s and nj, corresponding 
to the short and long period components of libration. In the 
ERTBP, however, the point L4 has an elliptic motion itself, 
similar to the elliptic motion of the primaries, and a corre- 
sponding mean motion n. This frequency combines with the 
frequencies of the short and long period librations, and thus 
motions around L4 have four frequencies (and linear combi- 
nations of them) which for e — > approach the limit values 
n s , n — m, m, and n — n s . Taking n as unit fr equency, the 
norm alized frequencies n a and n\ are given by (|Rabelll970l . 
[19Z3): 



-(l + vT^t), 



2 

77-1 



where 
27 



4 + 



(1) 



(2) 



The frequencies n s and n\ are real until /1 ^ 0.04006 . . ., 
and when the equality is held, n s — n\ — l/y/2 is a double 
root. (Note that values of /x of this section can be simply 
transformed to the system, where the primaries have masses 
1 — fj, and fi, dividing them by 1 + p. This results in the 
familiar 0.03852 . . . value for the critical mass parameter.) 

The normalized frequencies n$, 1 — n\, n\, and 1 — n s 
are given as the function of the mass parameter fi in Fig. 
6. The frequency curves n$([i) and ni(/x) meet at the point 
B, the curves 1 — m and 1 — n s at the pont C, both for 
/j, = 0.04006, and the common frequencies are l/\/2 and 
1 — l/\/2, respectively. The curves 1 — n\ and n\ cross at 
the point A for /j, — 0.02944 and the common frequency is 
1/2. The points A, B, and C correspond to 1:1 resonances 
between the relevant frequencies. 
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Figure 6. Limit frequencies of librational motions in the elliptic 
restricted three-body problem for e — > 0. 

In the ERTBP, the four frequencies of librations a round 
Lj ca n be determined from an equation derived by iRabel 
l| 1973ft . which relates the frequency z, the modified mass 
parameter k (see Equation 2) and the eccentricity e: 

D(z, k, e) = A(z, k) + e 2 B(z, k) + e 4 C(z, k) = 0. (3) 

The functions A(z,k), B(z,k), and C(z,k) are given in the 
Appendix. Equation (3) is 20th order in z, 5th order in n, 
and 4th order in e. Still, this is an approximate equation, 
obtained from an infinite determinant of an infinite number 
of homogeneous linear equations for the coefficients of oscil- 
lating solutions around L4. Nevertheless, Equation (3) can 
be used to determine the frequencies for any given k (corre- 
sponding to_0<u^ 0.05) and for small and moderate val- 
ues of e. IRabel (|l973T l used Equation (3) for a fourth-order 
stability analysis of L4 in the ERTBP, and to determine the 
transition curve between stability and instability in the k, e 
plane. 

In general, Equation (3) has four different roots for z 
between and 1, z = 1 corres ponding to the unit frequency 
of the primaries. IRabel I] 1973ft determined two solutions for 
Equation (3): 

(i) e = 0, k = 3/16 {n = 0.02944), 2 = 1/2; 

(ii) e = 0, k = 1/4 (fi = 0.04006), z = 1/^2. 

The solution (i) corresponds to the point A in Fig. 6, while 
the solution (ii) to the point B. These are critical points, 
where two frequencies are equal, and in the k, e plane bound- 
ary curves between stability and instability are defined by 
these two fi xed frequencies. Starting from these solutions, 
iRabd (|l973t ) derived approximate analytical solutions from 
Equation (3) for the boundary curves in the form of series 
expansions, valid for moderate values of e. 

For z = 1/2, Rabe's approximate solution is of third- 
order: 

o /oo 

n = — ±- — e + 0.04165e 2 =p0.14926e 3 . (4) 
16 16 

Adding to this a numerically determined fourth-order term, 
0.07e 4 , the solution (4) is shown in Fig. 5(b) as the curve 
labelled by All. As can be seen, it gives quite well the V- 
shaped boundary curve between the stable and unstable re- 
gion of Z/4 for moderate eccentricities (e < 0.3). Along this 
curve, 1 — n\ = n\ = 1/2, thus the boundary curve is the 
place where there is a 1:1 resonance between the frequencies 



Table 1. Types of resonances 



type 


resonance 


A 


(1 - m) : ni 


B 


n s : m 


C 


(1-m) : (l-n B ) 


D 


n s : (1 - n\) 


E 


n s : (1 - n s ) 


F 


n\ : (1 - n s ) 



1 — ri\ and n\, and moreover a 2:1 resonance between the 
frequencies 1 (the frequency of the primaries) and n\. 

For z = l/x/2. lRabel(|l973ft derived a fourth-order solu- 
tion. This solution gives the right side boundary between the 
stable and unstable region of L4 in the k, e plane. However, 
it seems that the fourth-or der term is erron eous in Rabe's 
solution. For this boundarv lTschaunerl (|l97ll ) derived an ex- 
act analytical solution. Using Tschauner's solution, Rabe's 
solution can be corrected, resulting in: 



This solution is shown in Fig. 5(b) as the curve Bll. Along 
this curve, n s = n\ = l/\/2, thus this boundary curve is the 
place where there is a 1:1 resonance between the frequencies 
n s and n\. As can be seen, the theoretical Bll curve fits 
very well the numerically obtained boundary, except for low 
values of the eccentricity (e < 0.05). We do not have an 
explanation for this yet. 

Apart from the above two solutions, no other solution 
was given to Equatio n (3). This m ight be due to the alge- 
braic difficulties which lRabd |l97fjft noted in connection with 
an earlier, simpler form of the equation. However, Equa- 
tion (3) offers the possibility to compute the frequencies for 
given mass parameters and eccentricities. We suspect that 
the structures in Fig. 5(a) are connected with resonances be- 
tween the frequencies. We have already seen that the bound- 
ary curves between the stable and unstable regions of L4 are 
places of 1:1 resonances. Considering the frequency curves in 
Fig. 6, several other types of resonances are possible. Table 
1 lists the possible pairs of the frequencies with a suggested 
name for the type of the resonance. 

Equation (4) is an approximate solution for the A 1:1 
resonance, while Equation (5) is for the B 1:1 resonance 
(that is why they are labelled as All and Bll in Fig. 5b). 
As can be seen from Fig. 6, there is also a C 1: 1 reso nance 
with frequency 1 — 1/ V2- It was shown bv lRabe! (|l970h that 
due to the structure of Equation (3), if z is a solution, so is 

1 — 2. Since Equation (5) is a solution for z — 1/V2, thus it 
is also a solution for 2 = 1 — 1/V2. Therefore, Equation (5) 
is also an approximate solution for the C 1:1 resonance. 

Numerically, a 1:1 resonance solution can be determined 
from Equation (3) by fixing the frequency at the resonant 
value and solving the equation for e with a given k, and 
repeating the procedure for other values of k. Thus an e(«) 
curve is obtained corresponding to the 1:1 resonance. Since 
Equation (3) is fourth-order in e, reducible to second order, 
the solution is easy. 

Solutions for other resonances can be obtained in the 
following way. If there are two frequencies in resonance, say 

2 and rz, where r is the ratio of the resonance, then both 
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Table 2. Initial points of the resonant curves at e = 



type 


resonance 


frequency 




ft 


A 


5:1 


0.83333333 





00403325 


A 


1:1 


0.80000000 





00575455 


A 


3:1 


0.75000000 





00883461 


B 


3:1 


0.94868330 





01370120 


A 


2:1 


0.66666667 





01507644 


D 


3:2 


0.92307692 





01940524 


B 


11:5 


0.91036648 





02195743 


B 


2:1 


0.89442719 





02489879 


A 


6:5 


0.54545455 





02554130 


A 


1:1 


0.50000000 





02943725 


B 


1:3 


0.80000000 





03668353 


B 


1:1 


0.70710678 





04006421 



frequencies have to satisfy Equation (3). Then we have two 
equations 

D(z, k, e) = 0, D(rz, k, e) = 0. (6) 

For a given value of k, we express the solution for e from 
both equations, containing z as parameter, and search for 
such value of z, for which the two solutions for e become 
equal. (In the formula for the solution of the second order 
equation giving e, the square root must be taken with nega- 
tive sign to obtain good solution.) Repeating the procedure 
for other values of k, the curve e(n) of the resonance r is ob- 
tained. Note that the resonant frequency is changing along 
this curve. 

Using the above method, we computed several resonant 
solutions. The corresponding resonant curves are shown in 
Fig. 5(b). Table 2 gives the initial points of these curves 
(mass parameter /j, at e = 0) and one of the initial resonant 
frequencies (dividing it by the resonant ratio, the other fre- 
quency can be obtained). Table 2 also contains the initial 
points of the solutions All, Bll given by Equations (4) and 
(5). Since Equation (3) is fourth-order in e, it gives satis- 
factory solution until about e < 0.3. Thus we computed the 
resonant curves below this limit. 

In Fig. 5(b), the two most prominent minimum zones 
are connected to two 2:1 resonances, A 2:1 and B 2:1. The 
B21 resonant curve goes sharply in the middle of the min- 
imum zone, which is bordered by the curve of the A 6:5 
resonance. The minimum zones on the left of the figure are 
connected mainly with high order A type resonances. 

Close to the A21 resonant curve, there goes the B31 
curve and they meet at the point /i c = 0.0103576, e c = 
0.265615. A third resonant curve, D32 goes exactly through 
this point, indicating multiple resonances. Noting, that the 
frequencies change along a resonant curve, it is convenient 
to denote by z a the frequency whose limit frequency is 1 — 
n\. Then the three frequencies at the crossing point of the 
three resonant curves are: n s = 0.9629994, m = 0.3209998, 
z s = 0.6419995. Indeed, n s /m = 3/1, z s /m = 2/1, and 
n s /zs = 3/2. There is also a Laplace-type relation between 
the three frequencies: z s + n\ — n s = 0. Computations show 
that the size of the stability region around L4 for the fj, c and 
e c value of the cross point is very reduced, it is less than 2° 
in synodic longitude. 

The starting points of the B31 and B21 resonant 
curves are at fi = 0.01370120 and /i = 0.02489879, cor- 




0.01 0.02 0.03 0.04 0.05 

Voss parameter 

Figure 7. Ratio of stable orbits depending on the mass parameter 
and eccentricity. 

responding to 0.013516016 and 0.024293897 in the usual 
mass parameter system. Th e se ar e the values for which 
iDeprit fc Deprit-Bartholomel (1 19671 ) showed the non-linear 
instability of L4 in the circular restriced three-body prob- 
lem. Interestingly, as can be seen from Fig. 5, along the A21 
resonant curve the size of the stability region around L4 
seems less extended than along the nearby B 3:1 resonance, 
though L4 is non-linearly stable at its e = starting point. 

On the right of the figure, the most prominent minimum 
zone is connected with the B 4:3 resonance. Interestingly, 
this zone begins at e « 0.1 and not at e = 0. 

It is worth to compare Fig. 5(a) with the results de- 
scribed in Section 2 concerning the maps of dynamical sta- 
bility. For each pair of the mass parameter and the gi- 
ant planet's eccentricity, we computed the LCI for about 
1.5 x 10 4 orbits by changing the semi-major axis and the 
orbital longitude of the Trojan planet. From this we can cal- 
culate the ratio of the number of stable orbits with respect 
to all the investigated orbits. This ratio can be assigned to 
each pair of fi and e, and in this way a map can be obtained 
showing the ratio of stable orbits. The results are shown in 
Fig. 7 (note that we extended the computations of Section 2 
up to /1 = 0.05). Though the resolution of this figure is much 
less than in Fig. 5(a) (and the mass parameters are slightly 
different), still the basic structure is similar. The minimum 
zones, identified in Fig. 5(b) as the A 2:1, A 3:1 and A 4:1 
resonances, appear also in Fig. 7. The zone of the A 2:1 res- 
onance at /i ~ 0.015 indicates an especially sharp minimum. 
Fig. 7 also comfirms the more probable occurance of stable 
orbits around L4 for small values of the mass parameter and 
the eccentric ity. We note that Fig . 7 is in good agreement 
with Fig. 6 of lErdi fc Sandorl (|2005h obtained by the method 
of the RLI. 



5 CONCLUSIONS 

In studying the possibility of Trojan exoplanets, the size of 
the stability region around the Lagrangian points L4 and 
I/5 is crucial. We determined the size distribution of the 
stability region around L4 in the elliptic restricted three- 
body problem depending on the mass parameter and the 
eccentricity. The size distribution has a fine structure with 
several minimum zones. We have shown that these minimum 
zones are connected with resonances between the frequencies 
of librational motions around L4. The situation is analogous 
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to the Kirkwood gaps at mean motion resonances in the 
distribution of the main belt asteroids of the Solar system. 

One main minimum zone, the region of the A 2:1 reso- 
nance is situated around /i = 0.015, which corresponds to a 
celestial body of about 15 Jupiter-mass around a solar-mass 
star. This could be called brown dwarf limit in the sense that 
it is not probable that brown dwarfs of this size can possess 
Trojan companions due to the reduded size of the stability 
region. However, the stable region above the brown dwarf 
limit in Fig. 5 suggests the possibility of Trojan planets ac- 
companied to larger objects. The regions of the higher order 
A 5:1, A 4:1, and A 3:1 resonances also indicate mass pa- 
rameter values (0.004, 0.006, 0.009) and correspondig eccen- 
tricities for which single planet exoplanetary systems hardly 
can have Trojan planets. Circumstances are most favorable 
for Trojan planets to exist, when both the mass parame- 
ter and the giant's planet eccentricity are small. However, 
higher order resonances may also affect this region. 
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APPENDIX A: COEFFICIENTS OF RABE'S FREQUENCY EQUATION 

A = 576z 4 - 28002 6 + 55962 s - 5925z 10 + 3581z 12 - 12502 14 + 246z 16 - 25z ls + z' 
+k (-576z 2 + 88O2 4 - 28z 6 - 295z 8 + 24z 1() + IO2 12 - 2O2 14 + 52 16 ) 
+k{12QQz 2 - I86O2 4 + 230z 6 + 3302 8 + 902 10 + IO2 12 ) 
+k 3 (144 - 202 2 + 5902 4 + 1402 6 + 102 s ) 
+k 4 (24 + 552 2 + 52 4 ) + k 5 , 



B = 138242 4 - 589922 s + 1029362 8 - 953192 10 + 5II682 12 - 163862 14 + 30602 16 - 3032 18 + 12 2 ; 
+k (-158402 2 + 147042 4 + 30582 6 + 11562 8 - 44562 10 + 18562 12 - 5622 14 + 842 16 ) 
+«: 2 (332082 2 - 276782 4 + 21242 6 + 42382 8 + 8882 10 + I8O2 12 ) 
+k 3 (3912 - 13702 2 + 71122 4 + 16742 6 + 1562 8 ) 
+k 4 (326 + 5272 2 +482 4 ), 



r 7iiofi,4 00.100.6 , 1249017 8 1188315 10 2551517 12 202743 14 75465 16 18 597 20 

C = 7116bz — 66b±22z H z z -\ z z H — z — 18742 H —z 

2 2 8 2 4 8 

( ,.o„™ 2 4 166205 6 272839 8 234923 10 „„„„„ 12 11251 14 5529 16 \ 

+k -1537202 2 + II88II2 4 H 2 6 H — 2 8 2 i0 + 229182 12 ^2 i4 + — — 2 16 

V 4 8 4 2 8 / 

W (31392 + 2S56732 2 - 2089062 4 + ™™f + ^2 8 + + 

^ ( 33952 + 55642 2 + + ^/ + 

V 8 4 8 / 

'1589 5645 2 249 4 \ 



+ 



